function [varargout]=Equ_NK(PP,SS,order)

%'Eta_Z','Eta_Zp','Eta_Z_N','Eta_Z_Np','Eta_Z_H','Eta_Z_Hp','Eta_M_dom','Eta_M_domp',...
 %'Eta_M_ext','Eta_M_extp','Eta_Y_H','Eta_Y_Hp','Eta_p_F','Eta_p_Fp','N_N','N_H',...
 %'w_N','w_H','mc_N','mc_H','Y_N','Y_H','Profit_N','Profit_H',...
 %'PAC_N','PAC_H','Y_Np','Y_Hp','K_N','K_H','K','Kp',...
 %'KI','KIp','Q_P','Q_Pp','R','Rp','L_N','L_H',...
 %'C_N','C_T','C_H','C_F','p_T','p_N','p_H','p_F',...
 %'Pir','Pirp','Pir_N','Pir_H','Pir_F','Lag_p_T','Lag_p_N','Lag_p_H',...
 %'Lag_p_F','Lag_p_Tp','Lag_p_Np','Lag_p_Hp','Lag_p_Fp','p_Np','p_Hp','Pir_Np',...
 %'Pir_Hp','Bp_dom','ir_dom','ir_ext','rr_dom','rr_ext','ir','TAU',...
 %'B','T','ir_s','Pir_H_s','Pir_F_s','dE','dEp','Y',...
 %'C','Profit','PAC','Lag_B','Lag_Bp','Lag_ir','Lag_irp','Lag_Eta_p_F',...
 %'Lag_Eta_p_Fp',
 
 PP_RHO_Z=PP.RHO_Z;
PP_RHO_Z_N=PP.RHO_Z_N;
PP_RHO_Z_H=PP.RHO_Z_H;
PP_RHO_M_dom=PP.RHO_M_dom;
PP_RHO_M_ext=PP.RHO_M_ext;
PP_RHO_Y_H=PP.RHO_Y_H;
PP_RHO_p_F=PP.RHO_p_F;
PP_KShare_N=PP.KShare_N;
PP_KShare_H=PP.KShare_H;
PP_Ret2Sca_N=PP.Ret2Sca_N;
PP_Ret2Sca_H=PP.Ret2Sca_H;
PP_DELTA=PP.DELTA;
PP_CapAdjCost=PP.CapAdjCost;
PP_EPSILON_N=PP.EPSILON_N;
PP_EPSILON_H=PP.EPSILON_H;
PP_THETA_N=PP.THETA_N;
PP_THETA_H=PP.THETA_H;
PP_Cons_FracT=PP.Cons_FracT;
PP_Cons_FracH=PP.Cons_FracH;
PP_Cons_ElasTN=PP.Cons_ElasTN;
PP_Cons_ElasHF=PP.Cons_ElasHF;
PP_Taylor_Pi=PP.Taylor_Pi;
PP_Taylor_Pi_H=PP.Taylor_Pi_H;
PP_Taylor_Pi_N=PP.Taylor_Pi_N;
PP_Taylor_dE=PP.Taylor_dE;
PP_Taylor_ir=PP.Taylor_ir;
PP_Flag_dE=PP.Flag_dE;
PP_Flag_Price=PP.Flag_Price;

 SS_Z_N=SS.Z_N;
SS_Z_H=SS.Z_H;
SS_N_N=SS.N_N;
SS_N_H=SS.N_H;
SS_w_N=SS.w_N;
SS_w_H=SS.w_H;
SS_mc_N=SS.mc_N;
SS_mc_H=SS.mc_H;
SS_Y_N=SS.Y_N;
SS_Y_H=SS.Y_H;
SS_Profit_N=SS.Profit_N;
SS_Profit_H=SS.Profit_H;
SS_PAC_N=SS.PAC_N;
SS_PAC_H=SS.PAC_H;
SS_C_T=SS.C_T;
SS_C_N=SS.C_N;
SS_C_H=SS.C_H;
SS_C_F=SS.C_F;
SS_L_N=SS.L_N;
SS_L_H=SS.L_H;
SS_p_T=SS.p_T;
SS_p_N=SS.p_N;
SS_p_H=SS.p_H;
SS_p_F=SS.p_F;
SS_Pir=SS.Pir;
SS_Pir_N=SS.Pir_N;
SS_Pir_H=SS.Pir_H;
SS_Pir_F=SS.Pir_F;
SS_Bp_dom=SS.Bp_dom;
SS_ir_dom=SS.ir_dom;
SS_ir_ext=SS.ir_ext;
SS_ir=SS.ir;
SS_TAU=SS.TAU;
SS_B=SS.B;
SS_T=SS.T;
SS_Y=SS.Y;
SS_C=SS.C;
SS_Profit=SS.Profit;
SS_PAC=SS.PAC;
SS_G=SS.G;
SS_C_H_s=SS.C_H_s;
SS_KI=SS.KI;
SS_K=SS.K;
SS_K_N=SS.K_N;
SS_K_H=SS.K_H;
SS_R=SS.R;
SS_Q_P=SS.Q_P;
SS_ir_s=SS.ir_s;

 Eta_Z=0;
Eta_Zp=0;
Eta_Z_N=0;
Eta_Z_Np=0;
Eta_Z_H=0;
Eta_Z_Hp=0;
Eta_M_dom=0;
Eta_M_domp=0;
Eta_M_ext=0;
Eta_M_extp=0;
Eta_Y_H=0;
Eta_Y_Hp=0;
Eta_p_F=0;
Eta_p_Fp=0;
N_N=0;
N_H=0;
w_N=0;
w_H=0;
mc_N=0;
mc_H=0;
Y_N=0;
Y_H=0;
Profit_N=0;
Profit_H=0;
PAC_N=0;
PAC_H=0;
Y_Np=0;
Y_Hp=0;
K_N=0;
K_H=0;
K=0;
Kp=0;
KI=0;
KIp=0;
Q_P=0;
Q_Pp=0;
R=0;
Rp=0;
L_N=0;
L_H=0;
C_N=0;
C_T=0;
C_H=0;
C_F=0;
p_T=0;
p_N=0;
p_H=0;
p_F=0;
Pir=0;
Pirp=0;
Pir_N=0;
Pir_H=0;
Pir_F=0;
Lag_p_T=0;
Lag_p_N=0;
Lag_p_H=0;
Lag_p_F=0;
Lag_p_Tp=0;
Lag_p_Np=0;
Lag_p_Hp=0;
Lag_p_Fp=0;
p_Np=0;
p_Hp=0;
Pir_Np=0;
Pir_Hp=0;
Bp_dom=0;
ir_dom=0;
ir_ext=0;
rr_dom=0;
rr_ext=0;
ir=0;
TAU=0;
B=0;
T=0;
ir_s=0;
Pir_H_s=0;
Pir_F_s=0;
dE=0;
dEp=0;
Y=0;
C=0;
Profit=0;
PAC=0;
Lag_B=0;
Lag_Bp=0;
Lag_ir=0;
Lag_irp=0;
Lag_Eta_p_F=0;
Lag_Eta_p_Fp=0;
if order==0
res=ones(59,1);
res(1)=Eta_Z*PP_RHO_Z - Eta_Zp;
res(2)=Eta_Z_N*PP_RHO_Z_N - Eta_Z_Np;
res(3)=Eta_Z_H*PP_RHO_Z_H - Eta_Z_Hp;
res(4)=Eta_M_dom*PP_RHO_M_dom - Eta_M_domp;
res(5)=Eta_M_ext*PP_RHO_M_ext - Eta_M_extp;
res(6)=Eta_Y_H*PP_RHO_Y_H - Eta_Y_Hp;
res(7)=Eta_p_F*PP_RHO_p_F - Eta_p_Fp;
res(8)=SS_Z_N*exp(Eta_Z + Eta_Z_N)*((SS_K_N*exp(K_N))^PP_KShare_N*(SS_N_N*exp(N_N))^(1 - PP_KShare_N))^PP_Ret2Sca_N - SS_Y_N*exp(Y_N);
res(9)=SS_Z_H*exp(Eta_Z + Eta_Z_H)*((SS_K_H*exp(K_H))^PP_KShare_H*(SS_N_H*exp(N_H))^(1 - PP_KShare_H))^PP_Ret2Sca_H - SS_Y_H*exp(Y_H);
res(10)=- SS_mc_N*exp(mc_N) - (SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(PP_Ret2Sca_N*SS_Y_N*(PP_KShare_N - 1));
res(11)=(SS_R*SS_K_N*exp(-Y_N)*exp(K_N + R))/(PP_KShare_N*PP_Ret2Sca_N*SS_Y_N) - SS_mc_N*exp(mc_N);
res(12)=- SS_mc_H*exp(mc_H) - (SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(PP_Ret2Sca_H*SS_Y_H*(PP_KShare_H - 1));
res(13)=(SS_R*SS_K_H*exp(-Y_H)*exp(K_H + R))/(PP_KShare_H*PP_Ret2Sca_H*SS_Y_H) - SS_mc_H*exp(mc_H);
res(14)=(Pir_Np*(SS_Pir_N + 1))/(SS_ir + 1) - Pir_N + (PP_EPSILON_N*(mc_N - p_N))/(PP_THETA_N*(SS_Pir_N + 1));
res(15)=((Pir_Hp - PP_Flag_Price*dEp*(SS_Pir_H + 1))*(SS_Pir_H + 1))/(SS_ir + 1) - Pir_H + PP_Flag_Price*dE*(SS_Pir_H + 1) + (PP_EPSILON_H*(mc_H - p_H))/(PP_THETA_H*(SS_Pir_H + 1));
res(16)=-PAC_N;
res(17)=-PAC_H;
res(18)=SS_Y_N*SS_p_N*exp(Y_N + p_N) - Profit_N - SS_PAC_N - SS_Profit_N - SS_R*SS_K_N*exp(K_N + R) - SS_N_N*SS_w_N*exp(N_N + w_N) - PAC_N;
res(19)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - Profit_H - SS_PAC_H - SS_Profit_H - SS_R*SS_K_H*exp(K_H + R) - SS_N_H*SS_w_H*exp(N_H + w_H) - PAC_H;
res(20)=(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN - SS_C_T*exp(C_T);
res(21)=- SS_C_N*exp(C_N) - (SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;
res(22)=(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF - SS_C_H*exp(C_H);
res(23)=- SS_C_F*exp(C_F) - (SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;
res(24)=PP_Cons_FracT*p_T - p_N*(PP_Cons_FracT - 1);
res(25)=PP_Cons_FracH*p_H - p_T - p_F*(PP_Cons_FracH - 1);
res(26)=exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1) - SS_Pir_N - Pir_N - 1;
res(27)=exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1) - SS_Pir_H - Pir_H - 1;
res(28)=exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1) - SS_Pir_F - Pir_F - 1;
res(29)=Eta_M_dom - ir + PP_Taylor_Pi*Pir + PP_Taylor_Pi_H*Pir_H + PP_Taylor_Pi_N*Pir_N + PP_Taylor_dE*dE + Lag_ir*PP_Taylor_ir;
res(30)=(B + SS_B)/(SS_ir + ir + 1) - SS_T - T - (Lag_B + SS_B)/(Pir + SS_Pir + 1) - SS_G + (SS_TAU + TAU)*(SS_L_H*SS_w_H*exp(L_H + w_H) + SS_L_N*SS_w_N*exp(L_N + w_N));
res(31)=-B;
res(32)=-T;
res(33)=ir - ir_dom;
res(34)=ir_s - ir_ext - dEp*(PP_Flag_dE - 1);
res(35)=ir_dom - Pirp - rr_dom;
res(36)=dEp - Pirp + ir_s - rr_ext;
res(37)=Eta_M_ext - ir_s;
res(38)=SS_C_H*exp(C_H) - SS_Y_H*exp(Y_H) + exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);
res(39)=Pir_F_s - Pir_F + dE;
res(40)=-Pir_H_s;
res(41)=Eta_p_F - Lag_Eta_p_F - Pir_F_s;
res(42)=(SS_KI*exp(KI - K))/SS_K - exp(Kp - K) - (PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K)^2)/2 - PP_DELTA + 1;
res(43)=SS_Q_P*exp(Q_P)*(PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K) + 1) - SS_p_F*exp(p_F);
res(44)=((SS_R*exp(Rp) - SS_Q_P*exp(Q_Pp)*(PP_DELTA + (PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K)^2)/2 + (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K - 1))*(Pirp + SS_Pir + 1))/(SS_ir + ir + 1) - SS_Q_P*exp(Q_P);
res(45)=SS_K_H*exp(K_H) - SS_K*exp(K) + SS_K_N*exp(K_N);
res(46)=B - Bp_dom + SS_B - SS_Bp_dom;
res(47)=SS_C_N*exp(C_N) - SS_Y_N*exp(Y_N);
res(48)=SS_L_N*exp(L_N) - SS_N_N*exp(N_N);
res(49)=SS_L_H*exp(L_H) - SS_N_H*exp(N_H);
res(50)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - SS_Y*exp(Y) + SS_Y_N*SS_p_N*exp(Y_N + p_N);
res(51)=PAC_H - PAC + PAC_N - SS_PAC + SS_PAC_H + SS_PAC_N;
res(52)=Profit_H - Profit + Profit_N - SS_Profit + SS_Profit_H + SS_Profit_N + SS_K*SS_R*exp(K + R) - SS_KI*SS_p_F*exp(KI + p_F);
res(53)=p_N - Lag_p_Np;
res(54)=p_T - Lag_p_Tp;
res(55)=p_H - Lag_p_Hp;
res(56)=p_F - Lag_p_Fp;
res(57)=B - Lag_Bp;
res(58)=ir - Lag_irp;
res(59)=Eta_p_F - Lag_Eta_p_Fp;
varargout{1}=res;
elseif order==1
fir_ord.Eta_Z=[PP_RHO_Z;0;0;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*((SS_K_N*exp(K_N))^PP_KShare_N*(SS_N_N*exp(N_N))^(1 - PP_KShare_N))^PP_Ret2Sca_N;SS_Z_H*exp(Eta_Z + Eta_Z_H)*((SS_K_H*exp(K_H))^PP_KShare_H*(SS_N_H*exp(N_H))^(1 - PP_KShare_H))^PP_Ret2Sca_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Zp=[-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_N=[0;PP_RHO_Z_N;0;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*((SS_K_N*exp(K_N))^PP_KShare_N*(SS_N_N*exp(N_N))^(1 - PP_KShare_N))^PP_Ret2Sca_N;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Np=[0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_H=[0;0;PP_RHO_Z_H;0;0;0;0;0;SS_Z_H*exp(Eta_Z + Eta_Z_H)*((SS_K_H*exp(K_H))^PP_KShare_H*(SS_N_H*exp(N_H))^(1 - PP_KShare_H))^PP_Ret2Sca_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Hp=[0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_dom=[0;0;0;PP_RHO_M_dom;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_domp=[0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_ext=[0;0;0;0;PP_RHO_M_ext;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_extp=[0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_H=[0;0;0;0;0;PP_RHO_Y_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_H*exp(PP_Cons_ElasHF*(p_F - p_H));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_Hp=[0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_p_F=[0;0;0;0;0;0;PP_RHO_p_F;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1];
fir_ord.Eta_p_Fp=[0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N_N=[0;0;0;0;0;0;0;-(PP_Ret2Sca_N*SS_N_N*SS_Z_N*exp(Eta_Z + Eta_Z_N)*exp(N_N)*((SS_K_N*exp(K_N))^PP_KShare_N*(SS_N_N*exp(N_N))^(1 - PP_KShare_N))^(PP_Ret2Sca_N - 1)*(SS_K_N*exp(K_N))^PP_KShare_N*(PP_KShare_N - 1))/(SS_N_N*exp(N_N))^PP_KShare_N;0;-(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(PP_Ret2Sca_N*SS_Y_N*(PP_KShare_N - 1));0;0;0;0;0;0;0;-SS_N_N*SS_w_N*exp(N_N + w_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_N_N*exp(N_N);0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N_H=[0;0;0;0;0;0;0;0;-(PP_Ret2Sca_H*SS_N_H*SS_Z_H*exp(Eta_Z + Eta_Z_H)*exp(N_H)*((SS_K_H*exp(K_H))^PP_KShare_H*(SS_N_H*exp(N_H))^(1 - PP_KShare_H))^(PP_Ret2Sca_H - 1)*(SS_K_H*exp(K_H))^PP_KShare_H*(PP_KShare_H - 1))/(SS_N_H*exp(N_H))^PP_KShare_H;0;0;-(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(PP_Ret2Sca_H*SS_Y_H*(PP_KShare_H - 1));0;0;0;0;0;0;-SS_N_H*SS_w_H*exp(N_H + w_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_N_H*exp(N_H);0;0;0;0;0;0;0;0;0;0];
fir_ord.w_N=[0;0;0;0;0;0;0;0;0;-(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(PP_Ret2Sca_N*SS_Y_N*(PP_KShare_N - 1));0;0;0;0;0;0;0;-SS_N_N*SS_w_N*exp(N_N + w_N);0;0;0;0;0;0;0;0;0;0;0;SS_L_N*SS_w_N*exp(L_N + w_N)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.w_H=[0;0;0;0;0;0;0;0;0;0;0;-(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(PP_Ret2Sca_H*SS_Y_H*(PP_KShare_H - 1));0;0;0;0;0;0;-SS_N_H*SS_w_H*exp(N_H + w_H);0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_N=[0;0;0;0;0;0;0;0;0;-SS_mc_N*exp(mc_N);-SS_mc_N*exp(mc_N);0;0;PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_H=[0;0;0;0;0;0;0;0;0;0;0;-SS_mc_H*exp(mc_H);-SS_mc_H*exp(mc_H);0;PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_N=[0;0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(PP_Ret2Sca_N*SS_Y_N*(PP_KShare_N - 1));-(SS_R*SS_K_N*exp(-Y_N)*exp(K_N + R))/(PP_KShare_N*PP_Ret2Sca_N*SS_Y_N);0;0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0];
fir_ord.Y_H=[0;0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;0;(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(PP_Ret2Sca_H*SS_Y_H*(PP_KShare_H - 1));-(SS_R*SS_K_H*exp(-Y_H)*exp(K_H + R))/(PP_KShare_H*PP_Ret2Sca_H*SS_Y_H);0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;0;0;0;0;0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0];
fir_ord.Profit_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.Profit_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.PAC_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.PAC_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.Y_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.K_N=[0;0;0;0;0;0;0;PP_KShare_N*PP_Ret2Sca_N*SS_K_N*SS_Z_N*exp(Eta_Z + Eta_Z_N)*exp(K_N)*((SS_K_N*exp(K_N))^PP_KShare_N*(SS_N_N*exp(N_N))^(1 - PP_KShare_N))^(PP_Ret2Sca_N - 1)*(SS_K_N*exp(K_N))^(PP_KShare_N - 1)*(SS_N_N*exp(N_N))^(1 - PP_KShare_N);0;0;(SS_R*SS_K_N*exp(-Y_N)*exp(K_N + R))/(PP_KShare_N*PP_Ret2Sca_N*SS_Y_N);0;0;0;0;0;0;-SS_R*SS_K_N*exp(K_N + R);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_K_N*exp(K_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.K_H=[0;0;0;0;0;0;0;0;PP_KShare_H*PP_Ret2Sca_H*SS_K_H*SS_Z_H*exp(Eta_Z + Eta_Z_H)*exp(K_H)*((SS_K_H*exp(K_H))^PP_KShare_H*(SS_N_H*exp(N_H))^(1 - PP_KShare_H))^(PP_Ret2Sca_H - 1)*(SS_K_H*exp(K_H))^(PP_KShare_H - 1)*(SS_N_H*exp(N_H))^(1 - PP_KShare_H);0;0;0;(SS_R*SS_K_H*exp(-Y_H)*exp(K_H + R))/(PP_KShare_H*PP_Ret2Sca_H*SS_Y_H);0;0;0;0;0;-SS_R*SS_K_H*exp(K_H + R);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_K_H*exp(K_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.K=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;exp(Kp - K) - (SS_KI*exp(KI - K))/SS_K - (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K;(PP_CapAdjCost*SS_KI*SS_Q_P*exp(Q_P)*exp(KI - K))/SS_K;-(PP_CapAdjCost*SS_KI^2*SS_Q_P*exp(2*KI - 2*K)*exp(Q_Pp)*(Pirp + SS_Pir + 1))/(SS_K^2*(SS_ir + ir + 1));-SS_K*exp(K);0;0;0;0;0;0;SS_K*SS_R*exp(K + R);0;0;0;0;0;0;0];
fir_ord.Kp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(Kp - K);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.KI=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(SS_KI*exp(KI - K))/SS_K + (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K;-(PP_CapAdjCost*SS_KI*SS_Q_P*exp(Q_P)*exp(KI - K))/SS_K;(PP_CapAdjCost*SS_KI^2*SS_Q_P*exp(2*KI - 2*K)*exp(Q_Pp)*(Pirp + SS_Pir + 1))/(SS_K^2*(SS_ir + ir + 1));0;0;0;0;0;0;0;-SS_KI*SS_p_F*exp(KI + p_F);0;0;0;0;0;0;0];
fir_ord.KIp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Q_P=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Q_P*exp(Q_P)*(PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K) + 1);-SS_Q_P*exp(Q_P);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Q_Pp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-(SS_Q_P*exp(Q_Pp)*(Pirp + SS_Pir + 1)*(PP_DELTA + (PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K)^2)/2 + (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K - 1))/(SS_ir + ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.R=[0;0;0;0;0;0;0;0;0;0;(SS_R*SS_K_N*exp(-Y_N)*exp(K_N + R))/(PP_KShare_N*PP_Ret2Sca_N*SS_Y_N);0;(SS_R*SS_K_H*exp(-Y_H)*exp(K_H + R))/(PP_KShare_H*PP_Ret2Sca_H*SS_Y_H);0;0;0;0;-SS_R*SS_K_N*exp(K_N + R);-SS_R*SS_K_H*exp(K_H + R);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_K*SS_R*exp(K + R);0;0;0;0;0;0;0];
fir_ord.Rp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(SS_R*exp(Rp)*(Pirp + SS_Pir + 1))/(SS_ir + ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.L_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_N*SS_w_N*exp(L_N + w_N)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_N*exp(L_N);0;0;0;0;0;0;0;0;0;0;0];
fir_ord.L_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*exp(L_H);0;0;0;0;0;0;0;0;0;0];
fir_ord.C_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_T*exp(C_T);0;(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF;-(SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_F*exp(C_F);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-(PP_Cons_FracT*PP_Cons_ElasTN*SS_C*SS_p_T*exp(C)*exp(p_T))/(SS_p_T*exp(p_T))^(PP_Cons_ElasTN + 1);0;(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));-(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));PP_Cons_FracT;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0];
fir_ord.p_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;(PP_Cons_ElasTN*SS_C*SS_p_N*exp(C)*exp(p_N)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^(PP_Cons_ElasTN + 1);0;0;1 - PP_Cons_FracT;0;exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;1;0;0;0;0;0;0];
fir_ord.p_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;-(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;0;PP_Cons_FracH;0;exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;-PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;0;0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;1;0;0;0;0];
fir_ord.p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;1 - PP_Cons_FracH;0;0;exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;-SS_p_F*exp(p_F);0;0;0;0;0;0;0;0;-SS_KI*SS_p_F*exp(KI + p_F);0;0;0;1;0;0;0];
fir_ord.Pir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;exp(p_N - Lag_p_N);exp(p_H - Lag_p_H);exp(p_F - Lag_p_F);PP_Taylor_Pi;(Lag_B + SS_B)/(Pir + SS_Pir + 1)^2;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pirp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-1;0;0;0;0;0;0;0;(SS_R*exp(Rp) - SS_Q_P*exp(Q_Pp)*(PP_DELTA + (PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K)^2)/2 + (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K - 1))/(SS_ir + ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;PP_Taylor_Pi_N;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;PP_Taylor_Pi_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_Tp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0];
fir_ord.Lag_p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0];
fir_ord.Lag_p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0];
fir_ord.Lag_p_Fp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0];
fir_ord.p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;(SS_Pir_N + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;(SS_Pir_H + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Bp_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-(B + SS_B)/(SS_ir + ir + 1)^2;0;0;1;0;0;0;0;0;0;0;0;0;0;-((SS_R*exp(Rp) - SS_Q_P*exp(Q_Pp)*(PP_DELTA + (PP_CapAdjCost*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K)^2)/2 + (PP_CapAdjCost*SS_KI*exp(KI - K)*(PP_DELTA - (SS_KI*exp(KI - K))/SS_K))/SS_K - 1))*(Pirp + SS_Pir + 1))/(SS_ir + ir + 1)^2;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0];
fir_ord.TAU=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H) + SS_L_N*SS_w_N*exp(L_N + w_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1/(SS_ir + ir + 1);-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;0;0];
fir_ord.T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;1;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dE=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Flag_Price*(SS_Pir_H + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_dE;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dEp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-(PP_Flag_Price*(SS_Pir_H + 1)^2)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1 - PP_Flag_dE;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y*exp(Y);0;0;0;0;0;0;0;0;0];
fir_ord.C=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN;-(SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Profit=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0];
fir_ord.PAC=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0];
fir_ord.Lag_B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1/(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_Bp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0];
fir_ord.Lag_ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_ir;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_irp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0];
fir_ord.Lag_Eta_p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_Eta_p_Fp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1];
varargout{1}=fir_ord;

 end
